Electrostatic Tuning of the Superconductor-Insulator Transition in Two Dimensions 



o 
o 



o 

O ; 
u • 

Oh' 

^3 ; 



X3 
O 

> 
O 

o 



I 

O 

o 



X 



Kevin A. Parendo, K. H. Sarwa B. Tan, A. Bhattacharya, M. Eblen-Zayas, N. E. Staley, and A. M. Goldman 
School of Physics and Astronomy, University of Minnesota, 116 Church St. SE, 
Minneapolis, MN 55455, USA 
(February 2, 2008) 



Superconductivity has been induced in insulating ultra- 
thin films of amorphous bismuth using the electric field ef- 
fect. The screening of the electron-electron interaction was 
found to increase with electron concentration in a manner cor- 
related with the tendency towards superconductivity. This 
does not preclude an increase in the density of states be- 
ing important in the development of superconductivity. The 
superconductor-insulator transition appears to belong to the 
universality class of the three dimensional XY model. 



Superconductor-insulator (SI) transitions in ultrathin 
films [1] are believed to be examples of quantum phase 
transitions (QPTs) [2]. Such transitions occur at zero 
temperature in response to the tuning of an external pa- 
rameter of a system that alters its ground state. Al- 
though many systems exhibit QPTs, including "^He ad- 
sorbed on random substrates, two-dimensional electron 
gases, and numerous complex strongly correlated elec- 
tron materials including high temperature superconduc- 
tors, the SI transition is one of the most fundamental be- 
cause of its apparent connection with the phase-number 
uncertainty relation of quantum mechanics [3] . The stan- 
dard approaches to tuning SI transitions involve increas- 
ing film thickness and applying perpendicular magnetic 
fields. Both of these can introduce complications. In- 
creasing film thickness changes carrier density but may 
also alter the disorder landscape. Applying perpendicu- 
lar magnetic fields introduces vortices and the complexity 
of vortex physics. 

Here we report the electrostatic tuning of the SI tran- 
sition in ultra-thin amorphous bismuth (a-Bi) films using 
the electric field effect. This tuning changes the carrier 
density without altering the disorder landscape. Mod- 
ulation of superconductivity using the field effect is an 
old idea [4], which has been revived by recent work on 
high temperature (cuprate) superconductors [5,6]. How- 
ever, these studies have not revealed significant features 
of the transitions between the various ground states. The 
difficulties of preparing high quality cuprate films with 
thicknesses on the order of the electrostatic screening 
length and chemically doped so that they are near a 
phase boundary has limited the usefulness of these stud- 
ies. These constraints are either irrelevant or easily sat- 
isfied in work on homogeneously disordered ultra-thin 
films. 

A single-crystal of (100) SrTiOs (STO) served as both 



the substrate and the gate insulator. Its unpolished sur- 
face was mechanically thinned [7], leaving parallel sur- 
faces approximately 50 /zm apart. Platinum was pre- 
deposited at 300 K forming a 1000 A thick gate electrode 
on the thinned back surface and 100 A thick measure- 
ment electrodes on the epi-polished front surface. The 
substrate was then placed in a dilution refrigerator/UHV 
deposition apparatus [8]. A 10 A thick underlayer of 
amorphous antimony (a-Sb) followed by subsequent lay- 
ers of a-Bi were deposited in situ under UHV conditions 
(~ 10~® torr) through shadow masks. The substrates 
were held at liquid-helium temperatures during deposi- 
tion. Films grown in this manner are believed to be ho- 
mogeneously disordered [9]. 

Resistances were measured by employing standard 
four-probe DC methods using a 1 nA current, a value 
well within the linear regime of the I-V characteristics. A 
voltage source provided the gate-film bias. All electrical 
leads were filtered at 300 K with RC filters to attenuate 
60 Hz noise and 7r-section filters to attenuate RF noise, 
and at the mixing chamber of the refrigerator with 2 m 
long Thermocoax cables [10] to attenuate GHz noise. The 
measurement and deposition processes were alternated, 
with the film always kept below 14 K in a UHV environ- 
ment, permitting the study of the evolution of electrical 
transport properties with thickness. 

In the sequence of films in which electrostatic gat- 
ing was studied, the film thickness was varied from 9.60 
A to 10.22 A in four steps. All films were insulating, 
with R(T) exhibiting a temperature dependence consis- 
tent with Mott variable range hopping. All curves exhib- 
ited resistance flattening for T < 60 mK. We will argue 
later that this behavior is not evidence of an intermedi- 
ate metallic regime. It is believed to be a consequence 
of the electrons in the film failing to cool despite signif- 
icant efforts to shield from external noise and thermal 
radiation. 

The 10.22 A thick film was the first in the sequence 
in which an electrostatically tuned SI transition (Fig. 1) 
could be induced. Resistance decreased with increasing 
positive gate voltage (Vg) at all temperatures, consistent 
with previous observations that the carriers in metallic a- 
Bi are electrons [11]. A logarithmic dependence on tem- 
perature became a better fit at about 7 V and continued 
to do so for all larger gate voltages at high temperatures 
despite the presence of superconducting behavior at low 
temperatures. 

One can qualitatively distinguish between insulating 
and superconducting ground states by examining the sign 
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of dR/dT at low temperatures. At low Vg, dR/dT < 0, 
suggesting an insulating state. With increasing Vg, 
dR/dT eventually changed sign, suggesting a supercon- 
ducting state. Zero resistance, within the limits of the 
scatter of the data, was observed for Vg > 38 V. The 
highest temperature at which zero resistance was ob- 
served increased monotonically with Vg, from 52 mK at 
38V to 59 mK at 42.5 V. At 42.5 V, all effects due to 
gating saturated, so R(T) for Vq > 42.5 V was identical 
to R(T) at 42.5V. 

The response of the resistance to Vq was stronger 
at low temperatures than at high temperatures, imply- 
ing that the high temperature conductance was affected 
much less by electrostatic charging than the supercon- 
ducting pairing mechanism. The major changes in R{T) 
in response to Vg appeared only as low temperature de- 
viations from the weakly insulating Vg = curve. Qual- 
itatively similar behavior has been found in the perpen- 
dicular magnetic field- tuned SI transition of In203 films 
[12], as well as in studies of the metal-insulator transition 
[13]. 

The inset of Fig. 1 shows the field effect for both signs 
of gate voltage at 65 mK for the 10.22 A thick film. Neg- 
ative Vg produces a small change in film resistance, while 
positive Vg causes a large decrease. At higher temper- 
atures, all R(T) curves for negative Vq largely coincide 
with the Vg = curve. Thus hopping transport is only 
changed a small amount. Similarly, for a 10.69 A thick 
film (not shown) with a 446 mK transition temperature, 
the asymmetry in the response for the two signs of gate 
voltage was evident in that for Vg = 50 V, T,. increased 
to 502 mK, while for Vq = —50 V, it decreased only to 
436 mK. This film was prepared in an earlier sequence of 
depositions in which the superconductor-insulator tran- 
sition was overshot. 

Measurements at low temperatures using thinned STO 
as a gate dielectric have demonstrated that electron 
transfer is approximately linear with gate voltage from 
zero at Vg = to approximately 3 x lO^"* cm^^ at Vq = 
50 V [7] and it does not saturate above that value. Thus, 
the observed saturation of film behavior cannot result 
from properties of the gate dielectric. The film's intrin- 
sic electron concentration is unknown, as the geometry 
in this particular set of studies was limited to parallel 
magnetic fields, precluding Hall effect measurements. 

The simplicity of electrostatic charging as a tuning 
parameter provides an opportunity to analyze data in 
ways not possible for either the thickness or magnetic 
field tuned SI transitions. The behavior of R(T) at high 
temperatures, as a function of gate voltage, can reveal 
features of the disordered 2D electron system relevant 
to its superconductivity. In 2D metals, in zero mag- 
netic field, the conductance is the sum of the classical 
Boltzmann conductance (Gs) and corrections due to co- 
herent backscattering/weak localization (Gwl) [14] and 
electron-electron interactions (Gee) [15]. Both produce 
a decrease in conductance that goes as the logarithm 
of temperature. This is based on a perturbative anal- 



ysis valid when kpl ^ 1, where kp is the 2D Fermi 
wavevector and I is the electronic mean free path. How- 
ever, it is possible that this analysis may work even 
when kpl ^ 1- (The value of kpl for these films, esti- 
mated from their high temperature conductivity is ap- 
proximately 1.6.) Thus, 

Gwl{T) + Gee(T) = [ap + (1 - \f*)\^MT) 

The weak localization correction is determined by the 
temperature exponent of the inelastic scattering time, 
p, and the localization parameter, a. The interaction 
term is controlled by a screening parameter, F* , which 
increases from to 0.866 as the electron-electron in- 
teraction changes from being unscreened to completely 
screened. The inset to Fig. 2 shows a plot of G vs. ln{T) 
for the Vg = 12.5 V curve of the 10.22 A thick film. This 
curve is nonlinear at low T because the system in tending 
towards superconductivity. For T >175 mK, the curve is 
linear in ln{T), and the coefficient of the logarithm de- 
pends on the values of ap and F* . Assuming the weak 
localization contribution, ap, to be constant for all Vq, 
by subtracting slopes of G{ln{T)) at nonzero Vg from 
the slope of G{ln(T)) at Vg = 0, one can determine the 
value of AF* = F*{Vg) - F*{Vg = 0). This function, 
which is then the increase in screening constant induced 
by gate voltage, relative to the screening constant intrin- 
sic to the ungated film, is shown in figure 2. F* increases 
monotonically by approximately 0.35 from Vg = to 
Vg = 42.5 V. F* will increase with this functional form, 
no matter the precise value of ap or F*{Vg = 0). (Note 
that in this analysis, we force a logarithmic fit to the 
exponential insulating curves. This is fairly accurate at 
temperatures higher than roughly 150 mK.) A positive 
magnetoresistance, logarithmic in B, that would allow 
determination of ap and of F*{Vg = 0) is expected in 
fields B » kuT/gnB- We observed positive magne- 
toresistance for all B at 300 mK at Vg = 10 V, but no 
specific field dependence could be determined. 

The increase in F* with Vg is correlated with the 
enhancement of the tendency towards superconductiv- 
ity. In conventional metals the latter is controlled by 
the competition between the Coulomb repulsion and 
the phonon-mediated attractive interaction. Increased 
screening leads to the suppression of Coulomb repul- 
sion, ultimately permitting superconductivity to develop. 
Once superconductivity is achieved, increased screening 
should move to higher temperatures. Saturation of 
both Tc and F* with Vg might not be unexpected in the 
limit of strong metallicity, as the density of states for a 2D 
metal is a constant independent of carrier concentration. 
The electron density determines the value of conductiv- 
ity at high temperatures. Since a major change in the 
low temperature behavior occurs without a substantial 
change in the resistance at higher temperatures, screen- 
ing, which changes significantly, and not the areal charge 
density, might be controlling the transition to supercon- 
ductivity. However, in the BCS model, the product of the 
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density of states and the interaction potential determines 
the transition temperature. Therefore, an increase in the 
density of states may also contribute to the inducing of 
superconductivity. 

The asymmetry of the response to the sign of Vq would 
suggest that the Fermi energy of the ungated film is close 
to the mobility edge. Increasing the electron concentra- 
tion with positive Vq would then increase the screening 
and the density of states significantly, whereas decreasing 
the electron concentration with negative Vq would result 
in small changes in the hopping transport, screening and 
the density of states. 

It is expected that conductivity data associated with 
continuous QPTs can be collapsed using a finite-size scal- 
ing analysis [2] . As data has been acquired at increasingly 
lower temperatures, such analyses have been observed to 
fail because of flattening of the curves of R(T). Scaling 
also appears to fail when curves of R(T) that exhibit the 
onset of superconducting behavior fan down from weakly 
insulating curves. In the present work, finite size scaling 
of R(T) data with Vq as the tuning parameter was suc- 
cessful because we confined the temperature range and 
used fine increments of values of the control parameter. 
All data for T < 60 mK was excluded, assuming that 
the observed fiattening of R(T) was due to a failure to 
cool the electrons. In the inset of Fig. 3, R{Vg) is plot- 
ted at fixed temperature for fifteen isotherms between 
65 mK and 100 mK. This yields a well-defined crossing 
point at Vg = H V, which is taken to be the critical 
value of the tuning parameter, Vgc ■ The critical resis- 
tance (Re) of 19,100 Q is large in comparison with values 
found for the SI transition with other tuning parameters. 
Plotting R/Rc versus \Vg - Vgc\T~^'"^ for 54 values 
of Vg between and 42.5 V, the product vz = 2/3 is 
found to yield the best collapse of data onto two separate 
branches, one for superconducting films and the other for 
insulating films. This is shown in Fig. 3. Assuming the 
value of z to be 1, the value ly = 2/3 is consistent with 
the universality classes of cither the three dimensional 
(3D) XY or inverse XY models. The latter is believed 
to describe the transition to superconductivity in 3D [16] 
and might be expected for a QPT in a system governed 
by the 2D XY model at nonzero temperature with z = 1 
[2]. This value of the exponent product is difii'erent from 
those found using thickness as a tuning parameter for a- 
Bi films [1] and magnetic field as a tuning parameter for 
InOx [17] and MoGe [18] films, but is consistent with a 
study of the magnetic field tuned transition in a-Bi films 
[19]. 

The breakdown of scaling when data above 0.1 K is in- 
cluded suggests that the critical regime only extends from 
zero to some low temperature on the order of 0.1 K. Ex- 
cept for the limited temperature range, the quality of the 
curve is as good as or better than that of earlier scaling 
analyses, from the highest value of R/Rc down to 0.6 and 
for a range of \Vg - Vgc\T~'^''"' extending from to over 
700. The analysis is shown down to R/Rc = to demon- 
strate how isotherms deviate from the superconducting 



scaling branch. The bottom curve is for the lowest tem- 
perature and appears to adhere to the scaling function. 
Curves generated from data taken at increased temper- 
atures deviate at higher values of R/Rc- This implies 
that scaling works best at the lowest temperatures. In 
this work, all R(T) curves tend toward becoming temper- 
ature independent below approximately 60 mK. Curves 
for which 11V< Vg<38V have nonzero resistance at 
the lowest measured temperature, while curves for which 
Vg > 38 V show zero resistance. The fact that all R(T) 
curves from both groups satisfy the same scaling function 
for the same vz suggests that these two regimes are not 
different. This is the main reason for our earlier assertion 
that in this instance the observed fiattening of resistance 
is not intrinsic. In other experimental situations, the 
flattening could be a real metallic phase. 

Any scaling analysis has a lower temperature bound. 
Thus, data that appear on the insulating branch may ac- 
tually belong on the superconducting one, as R(T) curves 
that appear to be insulating might actually be observed 
to be superconducting if lower temperatures became ac- 
cessible. The continuous nature of the scaling function 
at the zero of |Vg ^ VgcI would allow for a higher critical 
resistance than that identified here without change of the 
critical exponent product. 

The electrostrictive response of STO can be ruled out 
as an explanation of the observed effects because R{Vg) 
is asymmetric in Vg while the electrostrictive response of 
STO is symmetric [7]. 

The thinness (lOA) of the a-Sb layer rules out the pos- 
sibility that screening may be caused by an accumula- 
tion of a mobile charge at the interface between the STO 
substrate and the a-Sb layer, similar to the field-effect 
modulation found between STO and AI2O3 [21]. The 
fact that negative gate voltage causes little if any change 
in the film is different from what was found in previous 
work on a-Bi grown on a- Ge [22] . 

In summary, the electric field effect has been used to 
tune the 2D SI transition. The evolution of supercon- 
ductivity is correlated with an increase in screening. The 
size of the critical regime appears to be quite small. The 
exponent product vz ~ 2/3 is foimd, so the universality 
class of the transition appears to be that of the 3D XY 
model. This product is different from the value 1/2; ~ 4/3 
that has been found in numerous previous experiments, 
which has been suggested to be a signature of classical 
percolation in 2D [20] . It is possible that the underlying 
mechanism for this electrostatically tuned transition is 
different from those in studies that are tuned by chang- 
ing film thickness or applying magnetic fields. 
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FIG. 1. R(T) as a function of Vg for the 10.22 A film. From 
top to bottom, gate voltages are 0, 2.5, 5, 7, 8, 9.6, 11, 12.5, 
14.5, 17, 19.5, 24, 28, 33, 38, and 42.5 V. Approximately 40 
data curves have been omitted for clarity. Inset: Resistance 
versus gate voltage for the 10.22 A film at 65 mK. 



FIG. 2. Change in the screening parameter AF* versus Vg- 
Inset: Conductance versus ln(T) for the Vg = 12.5 V curve. 



FIG. 3. Scaling analysis of R{T) with Vg as tuning param- 
eter. The critical exponent product vz = 2/3. Inset: J?(Vg) 

for isotherms from 65 to 100 mK, yielding a distinct crossing 
point at the critical gate voltage of 11 V. 
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